Classification of dispersion equations for homogeneous dielectric-magnetic 

uniaxial materials 



Ricardo A. Depine^^'^'^), Marina E. Inchaussandague^^" '') and Akhlesh Lakhtakia^^^) 

GEA — Grupo de Electromagnetismo Aplicado, Departamento de Fisica, 
Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, 

Ciudad Universitaria, Pabellon I, 1428 Buenos Aires, Argentina 
CONICET — Consejo Nacional de Investigaciones Cientificas y Tecnicas, 
Rivadavia 1917, Buenos Aires, Argentina 
^'^^ CATMAS — Computational and Theoretical Materials Sciences Group, 
Department of Engineering Science and Mechanics, 
Pennsylvania State University, University Park, PA 16802-6812, USA 



ABSTRACT 

The geometric representation at a fixed frequency of the wavevector (or dispersion) surface uj{k) for 
lossless, homogeneous dielectric-magnetic uniaxial materials is explored, when the elements of the 
relative permittivity and permeability tensors of the material can have any sign. Electromagnetic 
plane waves propagating inside the material can exhibit dispersion surfaces in the form of ellipsoids 
of revolution, hyperboloids of one sheet, or hyperboloids of two sheets. Furthermore, depending 
on the relative orientation of the optic axis, the intersections of these surfaces with fixed planes of 
propagation can be circles, ellipses, hyperbolas, or straight lines. The obtained understanding is 
used to study the reflection and refraction of electromagnetic plane waves due to a planar interface 
with an isotropic medium. 

Key words: Anisotropy, Negative refraction. Elliptic dispersion equation. Hyperbolic dispersion 
equation. Uniaxial material, Indefinite constitutive tensor 

1 Introduction 

Recent developments in mesoscopic (i.e., structured but effectively homogeneous) materials have 
significantly broadened the range of available electromagnetic constitutive properties, thereby al- 
lowing the realization of solutions to Maxwell's equations which could have been previously regarded 
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as mere academic exercises. Materials having effectively negative real permittivity and permeability 
have been constructed from arrays of conducting wires ^4^ and arrays of split ring resonators 

Such composite materials — often called metamaterials — exhibit a negative index of refraction 
in certain frequency regimes Under these conditions, the phase velocity vector is in the opposite 
direction of the energy flux, for which reason they have been called negative-phase-velocity (NPV) 
materials j71IH]. 

NPV metamaterials synthesized thus far are actually anisotropic in nature, and any hypothesis 
about their isotropic behavior holds only under some restrictions on propagation direction and 
polarization state. In anisotropic NPV materials, the directions of power flow and phase velocity 
are not necessarily antiparallel but — more generally — have a negative projection of one on 
the other p|. Since the use of anisotropic NPV materials offers flexibility in design and ease of 
fabrication, attention has begun to be drawn to such materials jlUl llll [T21 113j . 

Natural crystals are characterized by permittivity and permeability tensors with the real part of 
all their elements positive, a fact that leads to dispersion equations in the form of closed surfaces. On 
the other hand, a relevant characteristic of NPV metamaterials is that the real parts of the elements 
of their permittivity and permeability tensors can have different signs in different frequency ranges. 
As an example, Parazzoli et al. [2] demonstrated negative refraction using s-polarized microwaves 
and samples for which the permittivity and permeability tensors have certain eigenvalues that 
are negative real. Under such circumstances, dispersion equations are topologically similar to 
open surfaces ^1]. Consequently, the intersection of a dispersion surface and a fixed plane of 
propagation may be a curve of an unusual shape, compared with its analogs for natural crystals. 
For example, extraordinary plane waves in a simple dielectric (nonmagnetic) uniaxial medium can 
exhibit dispersion curves which are hyperbolic, instead of the usual elliptic curves characteristic of 
natural uniaxial crystals jl4[ I15j . In recent studies on the characteristics of anisotropic materials 
with hyperbolic dispersion curves, new phenomenons have been identified, such as omnidirectional 
reflection — either from a single boundary jlUj or from multilayers jl6| — and the possibility of an 
infinite number of refraction channels due to a periodically corrugated surface jl7l I18j . 

In this paper, we are interested in studying the conditions under which the combination of 
permittivity and permeability tensors with the real parts of their elements of arbitrary sign, leads 
to closed or open dispersion surfaces for a homogeneous dielectric-magnetic uniaxial material. To 
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characterize this kind of material, four constitutive scalars are needed: 

• e|| and /iy, which are the respective elements of the relative permittivity and relative perme- 
ability tensors along the optic axis; and 

• e_L and which are the elements of the two tensors in the plane perpendicular to the optic 
axis. 

These scalars have positive real parts for natural crystals, but their real parts can have any sign 
for artificial (but still effectively homogeneous) materials. The dispersion equation for plane waves 
in such a material can be factorized into two terms, leading to the conclusion that the material 
supports the propagation of two different types of linearly polarized waves, called magnetic and 
electric modes HO] • 

The relative permittivity and permeability tensors, e and fi, are real symmetric when dissipation 
can be ignored. Then, each tensor can be classified as j21j : 

(i) positive definite, if all eigenvalues are positive; 

(ii) negative definite, if all eigenvalues are negative; and 

(iii) indefinite, if it has both negative and positive eigenvalues. 

Thus, the relative permittivity tensor is positive definite if e_L > and ey > 0; it is negative definite 
if < and ey < 0; and it is indefinite if e^ey < 0. In the present context, we exclude constitutive 
tensors with null eigenvalues. A similar classification applies to the relative permeability tensor. If 
both e and jl are positive definite, the material is of the positive-phase-velocity (PPV) kind. 

The plan of this paper is as follows. Considering the different possible combinations of e and fl, 
we show in Section |21 that magnetic and electric propagating modes can exhibit dispersion surfaces 
which are 

(a) ellipsoids of revolution, 

(b) hyperboloids of one sheet, or 

(c) hyperboloids of two sheets. 
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As a byproduct of our analysis, we also obtain different possible combinations of e and /2 that 
preclude the propagation of a mode — either electric, magnetic or both — inside the material. In 
Section 131 we study the intersection between the dispersion surfaces and a fixed plane of propagation 
that is arbitrarily oriented with respect to the optic axis. We show that, depending on the relative 
orientation of the optic axis, different dispersion curves, in the form of circles, ellipses, hyperbolas 
or even straight lines, can be obtained. Previous studies on dielectric-magnetic materials with 
indefinite constitutive tensors only considered planes of propagation coinciding with coordinate 
planes, thus failing to identify the singular case of linear dispersion equations. These results are used 
in Section 21 to discuss the reflection and refraction of electromagnetic plane waves due to a planar 
interface between a dielectric-magnetic uniaxial material and an isotropic medium. Concluding 
remarks are provided in Section |SJ An exp(— iwt) time-dependence is implicit, with oo as angular 
frequency, t as time, and i = a/— T. 

2 Dispersion surfaces 

The relative permeability and permittivity tensors of the anisotropic medium share the same optic 
axis denoted by the unit vector c, and their four eigenvalues are denoted by e_|_ y and In 
dyadic notation 



with / the identity dyadic. In this medium, two distinct plane waves can propagate in any given 
direction: 

(i) electric modes, with dispersion equation 




(1) 



k ■ e ■ k = fJ.±e±e 



(2) 



and 



(ii) magnetic modes, with dispersion equation 



k ■ jl ■ k = kQ fi±e±fi 



(3) 



Here k is the wavevector and ko denotes the free-space wavenumber. 



4 



We decompose the wavevector k = k\\+ k± into its components parallel (/cy ) and perpendicular 
{k±) to the optic axis. After taking into account that 

k-e-k = e±{kxcf + ei\{k-cf, (4) 

p)) for electric modes can be rewritten as 

k^ ku 

Analogously, © for magnetic modes can be expressed as 

k^ kl 

^± + -^ = kW. (6) 
Equations © and © have both the quadric form 

which displays symmetry of revolution about the fcy axis in three-dimensional /c-space. The pa- 
rameters A and B depend on the kind of mode (electric or magnetic) and their values determine 
the propagating or evanescent character of each mode and the geometric nature of the dispersion 
surface for propagating modes. 

One of the following conditions applies for a specific mode: 

(i) A > and i? > 0: the dispersion surface is an ellipsoid of revolution; 

(ii) ^ > and i? < 0: the dispersion surface is an hyperboloid of one sheet (Figure^); 

(iii) A <Q and B > the dispersion surface is an hyperboloid of two sheets (Figure^); 

(iv) ^ < and B < 0: the mode is evanescent. 

Depending on the particular combination of e and /i, we obtain from these conditions different 
dispersion surfaces. For example, the dispersion equations for electric and magnetic modes in 
natural crystals are both represented by eq. d with A > and i? > 0, a fact that leads to 
the known result that electric and magnetic modes have dispersion surfaces in the form of either 
prolate or oblate ellipsoids of revolution. The same result is obtained for metamaterials with both 
constitutive tensors negative definite. When the analysis is repeated for all possible combinations 
between the four constitutive scalars e^, ey, fi± and /Uy, the results summarized in Table 1 are 
obtained. 
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3 Intersection with a fixed plane of propagation 



In the previous section, by considering plane wave propagation in an unbounded medium, we found 
the various geometric forms of the dispersion surfaces. At a specularly flat interface between two 
half-spaces filled with linear homogeneous materials, the tangential components of the wavevectors 
of the incident, transmitted and reflected plane waves must all be equal, and consequently, they all 
must lie in the same plane that is orthogonal to the interface. This plane is the plane of propagation. 
Let us now investigate the kinds of dispersion curves obtained when dispersion surfaces of the kind 
identified in Section [51 intersect by a specific plane of propagation, arbitrarily oriented with respect 
to the optic axis c. 

Without loss of generality, let the xy plane be the fixed plane of propagation in a cartesian 
coordinate system; furthermore, let c = CxX + Cyy + CzZ and k = kxX + kyy. The dispersion equation 
0, for electric modes, can then be rewritten as the quadratic equation 



Mil kl + 2Mi2 kx ky + M22kl = F , 



(8) 



where 



Mil 

Mi2 

M22 



-L- (^.. 



e± + (eil - e±)cl 



(9) 



F = k^e\\e± fj,_L 

The dispersion equation (j^J for magnetic modes also has the same quadratic form, but now the 
coefficients Mn, M12, M22, and F are obtained by the interchange {ey < — > /iy, e± < — > fi±^ in Q. 



M 



(10) 



The symmetric matrix 

' Mil M12 
M12 M22 

corresponding to the quadratic equation Q is defined by its three elements. This matrix can be 
diagonalized by rotating the xy plane about the z axis by a certain angle, thereby eliminating the 
kxky term in (jSJ. With vi and V2 denoting the orthonormalized eigenvectors of the matrix M, we 
can write k = kivi + k2V2- Likewise, with 



Ai =e± + (e|| -e±)(c2 +4) 
A2 = e± 

denoting the eigenvalues of M, we get the dispersion curve 



(11) 



Ai kl + \2kl = F 



(12) 
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in the plane of propagation. 

The dispersion curves for the mode represented by (|12|) can be classified by analyzing the signs 
of Ai, A2 and F. In particular, 

(i) if Ai, A2 and F all have the same sign, then the dispersion curve in the fixed plane of 
propagation is an ellipse, with semiaxes along the directions vi and 1)2; 

(ii) if Ai and A2 have both the same sign, but F has the opposite sign, then the mode represented 
by ((T^ is of the evanescent kind; 

(iii) if Ai and A2 have opposite signs, then the dispersion curve is a hyperbola, with semiaxes 
along the directions vi and V2] 

(iv) if one eigenvalue is equal to zero and the other (nonzero) eigenvalue has the same sign as F, 
then the dispersion curve is a straight line, parallel to the eigenvector associated with the null 
eigenvalue. 

4 Illustrative numerical results and discussion 

To illustrate the different possibilities for the dispersion curves, let us present numerical results for 
the following two cases: 

Case I: e_\_ = —2.1, ey = 1.9, ^_\_ = 1.3 and /Uy = —1.6; 
Case II: = 2.1, ey = —1.9, //^ = —1.3 and /iy = 1.6. 

Both constitutive tensors thus are chosen to be indefinite. According to Table 1, the electric and 
magnetic modes for both Case I and Case II have dispersion surfaces in the form of one-sheet 
hyperboloids of revolution, whose intersections with fixed planes of propagation are circles, ellipses, 
hyperbolas or straight lines — depending on the orientation of c. 

Furthermore, to show the usefulness of our analysis in visualizing dispersion curves for boundary 
value problems, let us now consider that the anisotropic medium is illuminated by a plane wave 
from a vacuous half-space, the plane of incidence being the xy plane. In terms of (a) the angle Oc 



7 



between the optic axis and the y axis and (b) the angle '^c between the x axis and the projection 
of the optic axis onto the xy plane, the optic axis can be stated as 

c = X sin 9c cos fc + y cos 9c + z sin 9c sin ipc , (13) 

and the eigenvalues A|', corresponding to electric modes can be written as 

Af = e_L + (ey - e_L)(l - sin^ 9c sin^ ipc 

^2 



\E — c 



(14) 



For Case I, < 0, A|' = e_L < 0, whereas the sign of Xf depends on the optic axis orientation. 
From (|14() we conclude for the electric modes as follows: 



Af > if 



sin^ 9c sin^ ipc < ^- — , (15) 

and the dispersion curves are hyperbolas with semiaxes along the directions vf and ■0|'; 
Af = if 

sin^ 9c sin^ ipc = ^- — , (16) 

ey - e± 

and the dispersion curves are straight lines parallel to the direction associated with the eigen- 
vector vf; and 



Af < if 



sin^ 9c sin^ ipc > ^- — , (17) 

and the dispersion curves are ellipses with semiaxes along the directions of the eigenvectors 
vf and vf . 

The same conclusions hold for electric modes in Case II. 

Analogously, the eigenvalues A^, corresponding to magnetic modes are as follows: 

Af = /i± + (/xy - n±){l - sin^ 9c sin^ (fc) \ s^ 
- Hj_ J 

For Case I, F*^ > and A^ = /i_L > 0. From ^ we deduce that 
• Af < if 

sin^ 9c sin^ ipc < 

/^y - Ai± 



sin^ 9c sin^ tpc < — — , (19) 



and the dispersion curves are hyperbolas with semiaxes along the directions uf and {"f ; 



Af = if 



sin^ Of. s\v? LP(. 



Mil 



(20) 

Mil - M± 

and the dispersion curves are straight hnes parallel to the direction associated of the eigen- 



vector V 



1 1 



Af > if 



2 2 Mil 

sin 9c sin ipc > 



(21) 

Mil - M± 

and the dispersion curves are ellipses with semiaxes along the directions of the eigenvectors 
vf^ and . 

The same conclusions hold for magnetic modes in Case II. 

Let ipc > 0° so that the optic axis is not wholly contained in the plane of incidence. There 
exist critical values of 9c at which the dispersion curve change from hyperbolic/elliptic to ellip- 
tic/hyperbolic. By virtue of (|16|) . the critical value for electric modes is given by 

1/2 



sm ( 



(eil - e±) sin'' ipc 



(22) 



Likewise, the critical value 



sm 



Mil 



1/2 



(23) 



(M|| - M±) sin^ 

for magnetic modes emerges from (j20|) . Expressions (|22|) and (|23|) are valid for both Cases I and 
II. At a critical value of 9c, the dispersion curve for the corresponding mode is a straight line. 

Suppose ipc = 60°, so that 9^ = 52.73° and 9^^ = 59.06°. Then, for 9c = 9f the dispersion 
curves in the plane of incidence are straight lines (electric modes) and hyperbolas (magnetic modes); 
whereas for 9c = 9^^ , the dispersion curves are ellipses (electric modes) and straight lines (magnetic 
modes). 

In Figure [21 the reciprocal space maps for four different orientations of the optic axis are shown: 



9c = 20° (both dispersion curves hyperbolic), 

9c = 9f = 52.73° (electric type linear and magnetic type hyperbolic), 
9c = 55° (electric type elliptic and magnetic type hyperbolic), and 
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• 9c = = 59.06° (electric type elliptic and magnetic type linear). 

For 9c > 9^^ = 59.06°, modes of both electric and magnetic types have elliptic dispersion cm'ves — 
just as for a natm'al crystal (not shown). The light gray circle in Figure [21 represents the dispersion 
equation for plane waves in vacuum (the medium of incidence). 

For 9c = 20°, Figure indicates the nonexistence of real- valued ky in the refracting anisotropic 
medium for either the electric or the magnetic modes, the specific being indicated by a dashed 
vertical line in the figure. This is true for both Cases I and II, for any angle of incidence (with respect 
to the y axis), and for any incident polarization state; hence, the chosen anisotropic medium behaves 
as an omnidirectional total reflector [Tn|. As the present-day construction of NPV metamaterials 
is such that the boundary is periodically stepped it is worth noting that the introduction of 
a periodic modulation along the surface would subvert the omnidirectional reflector effect, since 
a periodic modulation allows for the presence of spatial harmonics with tangential components 
of their wavevectors that can now satisfy the required matching condition. Gratings of this kind, 
contrary to what happens for all gratings made of conventional materials, have been recently shown 
to support an infinite number of refracted channels |17l I18| . 

When 9c = 9c = 52.73° the dispersion equation for refracted modes of the electric type is linear. 
It is posible to find two wavevectors with real-valued components that satisfy the phase-matching 
condition (the so-called Snell's law) at the interface, one belonging to the upper straight line and 
the other to the lower straight line in Figure l^b. As the direction of the time-averaged Poynting 
vector associated with electric modes is given by j20j 

S = ^^{kx cfe-k , (24) 

we conclude that the refracted wavevectors on the upper straight line do not satisfy the radiation 
condition for Case I, whereas wavevectors on the lower straight line do not satisfy the radiation 
condition for Case II. 

The direction of S given by (|24j) for modes of the electric type is normal to the dispersion curves 
and points towards y < 0, as required by the radiation condition. Ray directions coincide with 
the direction of S. As for the parameters considered in our examples, the z component of the 
time-averaged Poynting vector does not vanish, the ray directions are not contained in the plane 
of incidence. The projections of the refracted rays onto the xy plane (indicated by little arrows in 
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the figures) are perpendicular to the straight hues and independent of the angle of incidence. 

For refracted modes of the magnetic type and for the angle of incidence (= sin~^ kx/ko) shown 
in Figure it is also posible to find two refracted wavevectors with real-valued components 
satisfying the phase-matching condition at the interface, one belonging to the upper hyperbola 
(not shown) and the other to the lower hyperbola. The time-averaged Poynting vector associated 
with the magnetic modes is given by 



Therefore, we conclude that wavevectors on the upper hyperbola do not satisfy the radiation condi- 
tion for Case II, whereas wavevectors on the lower hyperbola do not satisfy the radiation condition 



component in the z direction. Ray projections onto the xy plane (indicated by little arrows in the 
figures) are perpendicular to the hyperbolas. 

The interface for both Cases I and II acts as a positively refracting interface for modes of both 
types, in the sense that the refracted rays never emerge on the same side of the normal as the 
incident ray |22j . 

When the angle 6c is increased to 55° (Figure [2t)) the dispersion equation for the refracted 
modes of the magnetic type is still hyperbolic, but the dispersion equation for the electric type is 
elliptic. Again, for both electric and magnetic modes, is it possible to find two wavevectors with 
acceptable real-valued components. From ()24() . we conclude that refracted electric modes on the 
upper part of the ellipse correspond to Case II, whereas electric wavevectors on the lower part of 
the ellipse correspond to Case I. On the other hand, wavevectors for the refracted magnetic modes 
on the upper hyperbola do not satisfy the radiation condition for Case II, whereas wavevectors on 
the lower hyperbola do not satisfy the radiation condition for Case I, as can be deduced from ()25|). 

Ray projections onto the xy plane corresponding to the magnetic modes alone are shown in the 
figure, for the sake of clarity. For both Cases I and II and for refracted modes of the electric and 
magnetic types, the refracted rays never emerge on the same side of the y axis as the incident ray, 
just as for positively refracting interfaces. 

When Oc = Of = 59.06° (Figure Hi), the dispersion curves for the refracted modes of the 
electric type continue to be ellipses, but now the dispersion curves for the modes of the magnetic 




(25) 



for Case I. Ray directions coincide with the direction of S given bv l251 which again has a non-zero 
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type become straight lines. For the electric modes, the selection of the wavevectors is identical to 
that in Figure For the refracted magnetic modes, wavevectors on the upper straight line do not 
satisfy the radiation condition for Case II, whereas wavevectors on the lower straight line do not 
satisfy the radiation condition for Case I. 

Ray projections onto the xy plane for the refracted magnetic modes are also drawn in the figure. 
Again, for both Cases I and II the surface acts as a positively refracting interface for modes of both 
types. 

5 Concluding remarks 

This work focused on the geometric representation at a fixed frequency of the dispersion surface 
Lo{k) for lossless, homogeneous dielectric-magnetic uniaxial materials. To encompass both natural 
crystals and the artificial composites used to demonstrate negative refraction (metamaterials), we 
assumed that the elements of the permittivity and permeability tensors characterizing the material 
can have any sign. We showed that, depending on a particular combination of the elements of these 
tensors, the propagating electromagnetic modes supported by the material can exhibit dispersion 
surfaces in the form of (a) ellipsoids of revolution, (b) hyperboloids of one sheet, or (c) hyperboloids 
of two sheets. Intersections of these surfaces with fixed planes of propagation lead to circles, ellipses, 
hyperbolas or straight lines, depending on the relative orientation of the optic axis. This analysis 
was used to discuss the reflection and refraction of electromagnetic plane waves due to a planar 
interface with vacuum (or any linear, homogeneous, isotropic, dielectric-magnetic medium). 
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Table 1: Types of possible dispersion surfaces for different combinations between the eigenvalues 

e±) c^nd /i|| of the real symmetric tensors e and jl. The first symbol indicates the mode: E 

(electric) or M (magnetic). The second symbol indicates the geometrical form of the dispersion 

surface: e (ellipsoids of revolution), hi (hyperboloid of one sheet), /i2 (hyperboloid of two sheets). 
The symbol n indicates that the corresponding mode is of the evanescent (i.e., nonpropagating) 

kind. 
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Figure 2: Reciprocal space maps for Cases I and II, when = 60°. (a) = 20°, (b) Qc = 52.73°, 
(c) 6c = 55°, and (d) 6^, = 59.06°. The hght gray circle represents the dispersion equation for plane 
waves in the medium of incidence. 
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